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ABSTRACT 

A question that often arises about the principles under- 
lying the Brown- Twiss stellar interferometer is "How can 
interference phenomena occur with light produced by inco- 
herent sources?" Starting from the familiar interference 
pattern produced by two coherent sources, one can proceed 
in simple steps to a picture of two incoherent sources pro- 
ducing an interference pattern that moves about at random. 
This randomly moving pattern leaves behind a "footprint" 
in the form of the intensity autocorrelation function. This 
report describes how the autocorrelation function for an 
extended, incoherent source may be constructed. It is 
this function that is measured by the Brown-Twiss stellar 
interferometer. 
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AN AUTOCORRELATION APPROACH TO THE 
BROWN-TWISS INTERFEROMETER 


by 

Frank C. Jones 
Goddard Space Flight Center 


INTRODUCTION 

The physical principles underlying interference phenomena, such as the pattern produced by 
a diffraction grating or a double-slit interferometer, are understood by practically everyone who 
has studied physical optics. However, 10 or more years after the basic papers of Brown and Twiss 
(References 1 and 2) describing the principles and operation of their post-detection stellar inter- 
ferometers, much confusion still exists in the minds of many physicists concerning the fundamental 
concepts employed. This condition persists in spite of the considerable number of published papers 
(Reference 3) that have covered the phenomenon in great detail. 

There appear to be two main areas of confusion. The first concerns the correlation of arrival 
times of two separate photons at two separate detectors. This situation appears to contradict a 
basic principle of quantum electrodynamics: as stated by Dirac, (Reference 4) "Each photon then 
interferes only with itself. Interference between two different photons never occurs." 

This question has been answered quite clearly by Purcell (Reference 5) and Mandel (Reference 
6) to the effect that, while two photons cannot interfere with each other, they may be correlated by 
virtue of their being bosons and their tendency to clump together in phase space. It would be diffi- 
cult to improve on the clarity of the treatment of the quantum problem found in Reference 5 and 6. 
Therefore, only a brief discussion of this matter is given in Appendix A for background. In fact, 
this treatment leans heavily on that of Mandel (Reference 6). The primary result of this discussion 
is that the operation of the stellar interferometer can be understood entirely from the point of view 
of classical electromagnetic theory— a point that has been emphasized by Brown and Twiss (Refer- 
ence 2). In fact, the validity of the semiclassical approach to a wide variety of electromagnetic 
phenomena has been demonstrated recently by Mandel and Wolf (Reference 7). 

The second question is purely classical in nature: r, How is it possible to have any kind of 
interference phenomena associated with light from incoherent sources? ?t Although this question 
has been discussed extensively, it still remains somewhat of a puzzle to those not familiar with the 
literature on the coherence properties of radiation from incoherent sources. A thorough review of 
this subject and extensive references have been given by L. Mandel and E. Wolf (Reference 8). It 
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is to this question that the present paper is addressed. The following paragraphs describe the 
spatial dependence of intensity correlations, starting with the simple, first-order interference 
pattern produced by two coherent sources and proceeding the physical principles of the post- 
detection stellar interferometer of Brown and Twiss. 

SPATIAL DEPENDENCE OF INTENSITY FLUCTUATION CORRELATION 

The essential element in the post-detection interferometer is an electromagnetic radiation 
detector whose response is proportional to the intensity I of the radiation. This can be a radio 
receiver with a "square-law" detector circuit or, in the optical region of the spectrum, a photo- 
electric cell. In Appendix A it is shown that, although the photoelectric cell responds to individual 
photons, this fact can be ignored for the purpose of studying the statistical correlation between the 
output of two separate detectors, and its output can be assumed to be directly proportional to the 
classical intensity. 

In reality, however, detectors do not have zero resolution time and also do not have an output 
0(t ) that is directly proportional to the instantaneous intensity l(t); the appropriate quantity will 
be u(t), where 


u(t ) a I I( t ) dt . (1) 

J t -T 

In this expression, T represents the resolution time of the detector. If the coherence time At of 
the radiation is defined as the characteristic time during which the intensity varies only slightly, 
and if T «At , we have 


u(t)<*TI(t) . (2) 

This assumption is quite reasonable for the radio-interferometer (Reference 1) but not in the optical 
case (Reference 2). The complications involved in the case where T > At in no way change the basic 
principles to be considered, but they do make the treatment more difficult. For this reason, we 
shall assume that the output of any detector under discussion is directly proportional to the instan- 
taneous intensity of radiation falling upon it. 

In the operation of the post-detection interferometer, two detectors, A and B, are separated by 
a distance s . Their outputs are passed through filters, which remove the steady component and 
allow only the fluctuations to pass. The filter output is proportional to I( t ) - (i), where the brackets 
( ) represent a time average. The outputs from the two filters are multiplied together and time- 
averaged in a circuit called a correlator; the output of the correlator C AB is therefore 

Cab* - <I A <0>] [W <I B (0>J> 

= <Ia x b) - Oa> <Ib> ■ (3) 
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Figure 1 — Post- defection interferometer block diagram* 

Figure 1 is a schematic representation of this 
arrangement. 

It is desirable now to examine the rela- 
tionship between the correlation of intensity 
fluctuations in two separated detectors and the 
structure of the radiation source. We start with 
the simplest possible picture of two sources, 1 and 2, located at r l 
a monochromatic wave of frequency x and amplitudes k x and A 2 . At first, consider their ampli- 
tudes to be constant. However, each source will have a phase, t t ( t) and f 2 (t) , that will change 
with time in a random way, but with dr dt as 1 -t « . The important point to note is that, although 

the two sources have no fixed relative phase, at any given time a certain instantaneous phase dif- 
ference, : x (t) - i 2 (t), does exist. Neglect of this point causes much of the confusion concerning 
the possible existence of interference phenomena for incoherent sources. 


Figure 2— Instantaneous interference pattern 
produced by two sources. 


and r 2 (Figure 2), each radiating 


At the origin of the coordinate system, the amplitude from both sources is given by 


A - A. exp i 


h 


(t 


o] 


+ A 


2 exp i r 2 


ct r 


; 2 (t -r 2 c)] 


(4) 


In this expression, the amplitudes A t and A 2 are the amplitudes at the origin from source 1 
and source 2, respectively, and r x and r 2 are measured in wavelengths c for simplicity. 
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The corresponding intensity will be given by 


1 (t) - A x 2 + A 2 2 + 2 Aj A 2 cos [rj - r 2 + <f> x (t - r x /c) - <f > 2 (t - r a /c)] . (5) 

If we move to a point at a distance x away from the origin, keeping t fixed, and assume that 
x << r i» r 2 (this is extremely valid for earth observations of astronomical objects), the distances 
from the sources to the new position are given by 

r 1 ^ r 1 “ * x * r 2 ' % r 2 - ? 2 * x , 


where f = r/r , a. unit vector. 

The intensity at the point x (since we may safely take A x (x) = \ x etc.) is 

I (x, t) = Aj 2 + A 2 2 + 2 Aj A 2 cos^rj - r 2 - x * e 12 + $ 12 (x, t ) J , (6) 

where e 12 = r 1 - r 2 and, for small values, is just the angular separation between sources 1 and 2; 
$ 12 (x, t) is the new relative phase angle given by 

$ 12 (x, t) = 4> x (t - r j/c + r , • x/c) - 4> 2 (t “ r 2 /c + r 2 * x/c) . 

We may now ask to what extent Equation 6 represents a well defined interference pattern. Note 
first of all that the overall phase of the cosine function is unknown, since r x - r 2 (to within a 
wavelength) and the phase factor $ 12 are not known. However, the behavior of Equation 6 as a 
function of x is well defined if the phase angles do not change appreciably (if they do, they change 
in a random way). To assure this, the condition r, ■ x^? 2 * x « cAt must exist since 4>( t) is 
essentially constant for times that are small compared to At. This can be achieved by requiring 
f 1 • x = ? 2 * x = o; however, this would make x * e 12 = o, and we would have no dependence on x. 
The best compromise, therefore, is to set (x/2)*(r 1 + ? 2 ) = o, or in other words to stay on a plane 
perpendicular to the average direction of the sources (see Figure 2). In practice (References 1 
and 2) this is achieved by inserting a time delay in the circuitry of the detector at x equal to 
-St = (1/2) (r 1 + r 2 ) * x/c, which has the same effect on the arguments of 4> x and <£ 2 . 

Using either approach, the arguments of < j> x and 4> 2 niay be written as 

t - r 1 /c + ? 1 • x/c - (1/2) (r 1 + r 2 ) * x/c = t-r 1 /c+© 12 * x/2c 


and 


t - r 2 /c - e 12 * x/2c 
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respectively. We now require that e 12 • x « 2cAt. But we have required d<£/dt^ 1/At«<u, and, 
since co = c in our units, 2cAt » 1. Therefore, e 12 • x may vary through many times 2tt without 
appreciably changing <f> x or 4> 2 . 

With these considerations Equation 6 can be written 

1 (x, t) = Ij + I a + 2 yij I 2 cos^e 12 ' X + o 12 (t)J , (7) 

where we have dropped the constant (and unknown) factor r x - r 2 and have written l x and i 2 for 
A 2 and A 2 2 ? respectively. Thus we have a standard coherent interference pattern— but one that is not 
stationary. It jumps from place to place, staying in one spot only for times less than At. It 
would be extremely difficult to determine its structure ( and hence e 12 , which would be of astro- 
nomical interest) by taking its photograph— although this has been done for the interference fringes 
produced by the light from two lasers (Reference 9)— or by scanning it with a single phototube. On 
the other hand, it is important to realize that the spatial pattern represented by Equation 7 has a 
certain internal structure that is independent of its overall position. The internal structure of a 
function can be described in terms of its autocorrelation function (Reference 10), which is a meas- 
ure of how the values of a function at two points separated by a fixed amount are related. 

For a given function F(x) the autocorrelation function </>(S) is defined by 

/* x / 2 

lA(S) = JiSir F(x')F(x' + 8) dx' . (8) 

J -x/2 


It is easy to see that if F(x) is a periodic function with a period x, 


vKS) 


■+r 

•' X 


F(x')F(x' +S) dx # 


(9) 


which is independent of x. Forming this quantity for the intensity pattern (Equation 7) gives 


K 5) 


fX + X 

I(x',t)I(x'+6,t)dx' = (ij + I 2 ) 2 + 2 1 1 I 2 cos(e 12 ■ 5) , 


( 10 ) 


where 


X = 2t76 12 /\ 0 12 | 2 , 

Equation 10 is independent of x and the phase $ 12 (t). This quantity, viewed as a function of 8 , 
allows one to extract the physically interesting parameter e 12 , even though there is no 
information about the position of the interference pattern. 
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From our point of view, the significance of the quantity 0(6) lies in the fact that changing the 
phase angle o 12 is equivalent to changing (e 12 • x). Therefore, the integral over x in Equation 10 
can be replaced by an integral over o 12 : Thus, 

i r 277 

'/'(S) = I I(x', t) I(x' + 6, t) d* 12 . (11) 

•'0 

If we take a time-average of the integrand i(x, t) I (x + 6, t), the phase angle <P 12 will vary 
randomly and uniformly over the range 0 to 2?r, and we will have performed a Monte Carlo integra- 
tion of Equation 11. Therefore, 

<I(x, t)I(x+ S, t)> = '/'(S ) = (ij + I 2 ) 2 + 2 Ij i 2 cos (e 12 • s) . (12) 

Equation 12 demonstrates the underlying principle of the stellar interferometer of Brown and 
Twiss. It shows that two sources of radiation that have no fixed phase relationship can cause the 
intensity at two different points to be correlated. Furthermore, the correlation depends on the 
separation of the detectors 6 and the separation of the sources e 12 in a way that is very reminiscent 
of the interference pattern that would have been obtained had the sources been coherent. This is not 
surprising since the correlation pattern is closely related to the underlying interference pattern 
that is moving about too rapidly to be seen. The picture of the sources considered so far is very 
simple and not very realistic. In the following paragraphs, the picture will be refined to resemble 
a real astronomical object; the basic principle will remain the one expressed in Equation 12. 

REFINEMENTS OF THE MODEL 

The foregoing paragraphs described sources that radiated a scalar amplitude A. The descrip- 
tion therefore corresponded to a completely polarized source or to measurements made with aligned 
polarizing filters placed over the detectors. The treatment of partially polarized sources is rather 
complicated and will not be attempted here. However, the extension to completely unpolarized 
sources is almost trivial if one assumes that waves of opposite polarization are completely uncor- 
related and also remembers that they do not interfere with each other. Thus, for completely un- 
polarized sources, i H = i i = 1/2, and 

I(x, t) = Ij + l 2 + yij l 2 |cos [e 12 • X + ® 12( ,|, (t)] + cos [e 12 • X + <t> 12(i) ]| . (7') 

From this, 

<I(x, t) I(x + 8, t)> = (l x + I 2 ) 2 + I t I 2 cos (e l2 • 6) . (12') 
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Combining Equations 12 and 12 T into one equation gives 


<I(X, t) I(x + 8, t)) = (i x + i 2 ) 2 + rij i 2 cos (e 12 • s) , (12”) 

where y - 1 for completely unpolarized light, and y = 2 for completely polarized light. 

It is a straightforward matter to generalize Equation 12” for the case of n sources with direc- 
tion vectors f . ; thus 


<I(x, t) I(x + s, t)> 




N 


Z] C 1 ! «« (*ij 



(13) 


where 



r. • 6 

i 


= 0 


At this point, a rather peculiar fact becomes apparent. Until now the point oscillators have been 
described as radiating a constant-amplitude wave with only the phase varying at random. There- 
fore for a single source, i. = const., and (i 2 ) = (I.) 2 . However, allowing N of the point sources 
in Equation 13 to coincide at a point to form a single source (setting e .. = o), gives 


<I(x) I(x + 8)> 


<I 2 > 

it 4 




<i > 2 + t £ hh 


~1 


(14) 


This situation indicates that there is no longer a simple, point source with constant 
amplitude, but rather one whose intensity fluctuates in time so that (i 2 ) > (i) 2 . Realizing that 
any macroscopic radiation source (where macroscopic implies a dimension that is very large 
compared to the wavelength) is in fact composed of many microscopic oscillators, we see that 
our ”point” sources should be considered as being composed of many, small identical oscillators. 
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To answer the question of how many oscillators should contribute to one source, consider Equa- 
tion 14 for N equal oscillators where i. = (l)/N. Thus, 


(I 2 ) = (O 2 + T J2 li Ij 

i*j-l 

■ < 0 * [• i] 

= (I ) 1 (l *-5 - 3jsj) (15) 


Choosing N = 00 gives for any point oscillator the relationship 


Vi 2 ) = f 1 + t) (I )) 2 


(16) 


It may be stated that the picture of a point source composed of a very large number of 
oscillators, all oscillating with a steady amplitude but constantly changing phase, is not very 
realistic. This is a valid objection. However, in Appendix B, it is shown that a far more 
realistic model of a radiation source also yields Equation 16. This fact is not too surprising; 
it is merely one more example of the statistical law of large numbers. This law asserts that the 
probability distributions of a quantity that is a sum of N random quantities approaches a gaussian 
as N becomes infinite. In our case, this means that the statistical properties of the output from 
N oscillators become independent of the properties of the individual oscillators as N becomes 
very large. 

Thus in deriving the expression for the autocorrelation of the intensity (Equation 12) we should 
have taken time-averages of the individual terms, taking into account the fluctuations of i 1 and i 2 . 
Equation 12 M becomes 


<I(X, t) I(x + 6, t)> = <(l, + I 2 ) 2 ) + y <Ij I 2 ) cos(e 12 • s) 

and instead of Equation 13 we have 

N 

T 22 (li COs(9i i ' 6) 


(17) 



i^j =1 


21 < r i x i> ( x + ^ c ° s ( e ii • s )) 

\i=l / i*j=l 


(18) 
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Applying Equation 16 and noting that (I. I.) = <I.) (i.) , since the intensity fluctuations of two 

different sources are independent, gives 


(I(x, t). I(x + S, t)> = (l + 2) 

i = l 


+ £ f 1 + T cos ( e *i ' S )J ■ 


(19) 


Allowing these sources to coincide to form one source will give 

< i2> = { i + t) £ + ( 1 *%) £ (liKi,) 

i = 1 i^j = 1 

■ ( i4 ^) L ft>« = (“*)(£<■.>] 

i.j=l \i = l / 

= (i + t)«i» 2 • (20) 


which is the same as Equation 16 for a point source. This expression shows that an infinity of 
oscillators plus several other infinities of oscillators behaves like an infinity of oscillators. 

Equation 19 can be cast into a form that is quite suggestive; recalling that 0 .. = f . - ?., 


N 

<!(*) I(* + «)> = 22 ^ t 1 + ^ COS ( ? i ‘ 8 ■ *J ■ s 3 


I IN \ IN 

1 22 ^ ii >) + 2 " 21 ^ ^ exp ■ 6 ■ i? i ■ s ) 


i. j=l 


IN 

!>*> exp (ir i * s) 


= <i> 2 + $ 


2 


(21) 



The generalization from N discrete sources to a continuous distribution of sources is im- 
mediate. With the substitutions (i .) - I(f)dft, 2 - Equation 21 becomes 


<I(x) I(x + S)> 


<I> : 


y 

+ T 




F(f ) exp (if ■ 5 ) dH 


( 22 ) 


where F(f ) is the source brightness per solid angle in the direction f . 

Returning to Equation 3, we may identify l A with I(x), and l B with I(x+s). It is easy to 
verify that <l(x))= (i(x + 5))= (i); therefore Equation 22 can be substituted into Equation 3 
to obtain for the output of the correlator 


C AB 00 [<!(*) K* + 8)> - <!(*)> (l(* + 8)>] = T J 


F(r) exp (if ■ 6 ) dH 


(23) 


The quantity F(s) = f F(f) exp (if • s)dfi is the quantity determined with the Michelson interferom- 
eter. By determining this quantity as a function of s , the complete brightness distribution over the 
face of a star F(f ) may be determined, using the inverse Fourier transform 


F(r ) 


( 277 ) 


-2 


I 


F(6 ) exp (- if ■ 8 ) 


(24) 


With the post-detection interferometer, only |F(8 )| 2 can be determined, and there is no informa- 
tion about the phase of F(s). We are, therefore, unable to reconstruct F(f ) uniquely without making ' 
additional assumptions. 

The foregoing result has been obtained under the condition that d^/dt be small compared to 
the frequency This is equivalent to assuming that the source is quasi-monochromatie, since 
any change in 4> introduces additional frequencies to the wave, and that the condition d^/dt « ^ is 
the same as kco « co; in other words, we are dealing with a narrow band of frequencies. If this 
condition is relaxed, Equation 7 is no longer valid since the phase difference $ 12 cannot be considered 
as independent of x. Because of this, Equation 12 should be replaced by 

<I(x, t) I(x + 8, t)> = (l t + I 2 ) 2 + 2I t l 2 c(0 12 • 8) . (25) 


where 


c ( 0 12 • s) = 2 (cos[e i2 ■ x + $ 12 (x, t)J cos [e i2 ■ x + 0 12 -8 + 1> 12 (x + S, t ) 


(26) 
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If we write $ 12 (x + s , t) = <i> 12 (x, t) + A $ 12 (x + s, t), we can expand the cosine functions to 
obtain 


C ( 0 12 ’ 5 ) = ^ cos (<9i2 ' 6 ) [cos 2 (o 12 ■ x) (cos 2 <I> 12 COS A$ 12 

- cos $ 12 sin $ 12 sin A$ 12 ) + sin 2 (fl 12 * x) (sin 2 $ 12 cos A$ 12 

+ sin $ 12 cos <b 12 sin A3> 12 ) - cos (0 12 • x) sin (0 12 * x)( cos 2 $ 12 sin A$ 12 

- sin 2 $ 12 sin £& 12 + 2 sin $ 12 cos $ 12 cos A<£ 12 )J 

- sin (0 12 * s) [cos 2 (0 12 * x)(cos 2 3> 12 sin A$ 12 


+ cos <I> 12 sin $ 12 cos A$ 12 ) + sin 2 (0 12 * x) (sin 2 $ 12 sin A<$ 12 

- sin <J> 12 cos $ 12 cos A$ l2 ) + cos (0 12 ■ x) sin (0 12 * x) (cos 2 <E> t 2 cos A$ 12 

- sin 2 ^> 12 cos A$ 12 - 2 sin 2 cos $ 12 sin A<£ 12 )J^ , (27) 

where the arguments of <l> 12 and ao 12 are understood. We shall now assume that $ 12 performs a 

random walk so that the distribution of A$ 12 does not depend on $ 12 . The primary justification for 
such an assumption is that any natural model of a light source would have this property. 

Recall that a uniform distribution of <i> 12 implies 


(cos 2 ® 12 ) 


(sin 2 S> 12 ) 


X 

2 


and 


thus Equation 27 becomes 


^sin $ 12 cos $ 12 ^ ~ 0 


C (°12 ' 5 ) ' cos { e i2 ’ 6 ) ( cos A ® 12 ) " sin(0 12 • 5) (sin A<& j 2 ^ 
- R e [exp ( ie !2 ‘ S ) ( eX P ( iA ^X 2)) ] - 


(28) 
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where R e means the real part, and (exp (iA$ 12 can depend only on o l2 • s (not on x) since we 
have assumed throughout the discussion that the statistical properties of the source do not depend 
on time and the arguments of ^ 2 are just retarded time. 

The development from this point on proceeds exactly as before with the simple substitution of 
the function c(fl 12 * s) for cos (o l2 • 5 ). Recall that A$ 12 = A£ 2 - thus Equation 21 becomes 


<I(x)I(x + 6>) - R e {^(l^) [ l + "2~ ex P (i'i * 8 - ir j * s) (exp (iA0. - iA0.)) ] j 


J- - j _1 


IN 

A> 2 '* L <■■> exp (irj • s) (exp (iA<£.)) x (i.^ exp (- ir } • s) (exp (- iA<£.)) 
i , j = l 


■ <■> 


i > 2 +■% 


IN 

( eX P( iA ^i)) exp (iiv * s) 


(29) 


Before proceeding, pause and consider the function (exp (iA^>. . Forming the autocorrelation 

function of the amplitude A. ( t ) gives 

H T ) 3 (A.*(t)A(t+T)>/|A| 2 


exp ( i cor 


) (exp [iA<£(r)]) 


exp ( icor) \p' (r) 


(30) 


The Fourier transform of u>( r ) is 


?(«') 


f 

1 4>( r ) exp (ioj't) dr 


4> f O' + co) 


(31) 


where 4>' O' ) is the transform of ifj 1 (r). The function 4> (a; 1 ) is called the normalized power 
spectrum (Reference 10) of the amplitude A. (t) . It has this name because it represents the 
relative amount of power radiated in the frequency interval co 1 and w' + d co* (Reference 10). 
Therefore writing 


(exp [iA4>. (r)]^ 



\p* ( 00 ') exp (- ico* r) dco 1 , 
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I 


1 


and substituting this expression in Equation 29 gives 


<I(X) I(X + S)> - <1)2 



6 - ico* r i J \p/ (co* ) dco* 


(32) 


Abandoning the practice of writing lengths in units of the wavelength gives 5 - 6 A = Soj/c . We 
also have r. = r. • s/c and therefore can write 


<100 l(x + S)> - <l> 2 



(co - co* )/c] \jt/ (co 1 ) d co* 


(33) 


Assuming that all of the sources have the same spectral distribution, i.e., </>.' (co* ) = \p* (co* ) for 
all i , we may go to the limit of a continuously distributed source and write as the equivalent of 
Equation 23 


00 

/// 


F(r) exp [ir * 6 (co - co* )/c] \p* (to*) d co* dO 


00 

j jj F(r) exp [ir * 6o/ / c] ^(co* ) dco 1 di^ 


j : 


F(5oj* /c) tpito' ) d co * 


(34) 


This result is valid for any spectral distribution, assuming that the detector itself is not frequency- 
dependent. In any real situation, (co* ) should be multiplied by the frequency-response function of 
the detector since any real detector is equivalent to an ideal detector with a frequency filter placed 
in front of it. 

From Equation 34 it can be seen that if ^ (co* ) is sharply peaked about a particular frequency 
(quasi-monochromatic), we obtain the previous result, Equation 23. The finite bandwidth of the 
input radiation causes a loss of information in the sense that the function F (Bco/c) gets "smeared 
out" over the frequency distribution. For example, if the characteristic size of the source is 
0 , F will have structure over values of 8 of order c/co 0 or larger. If the spread in frequency 
A co is such that A (c/co 0) = (c/co 0) A co/co ^ c/co 0 or A co/co ^ 1 , the structure will be lost and 
even the overall size of the source cannot be determined with great accuracy. Therefore, even 
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though the initial requirement d<£/dt « a> is not strictly required, it certainly indicates the best 
operating conditions, and, if it is violated such that d<£/dt * it is very difficult to obtain any 
information about the source. 

CONCLUSIONS 

Beginning with the simple notion of an interference pattern that moves randomly about, one can 
proceed through a series of simple steps to the spatial correlation pattern of light from an extended 
source. This correlation pattern exists because so-called incoherent sources are not totally in- 
coherent. There are short periods of time during which some phase -relations exist between the 
various parts of the source. During this same short period, an interference pattern exists at the 
point of observation. Since these phase-relations between the various parts of the source do not 
persist, this interference pattern is constantly shifting and changing. However, during this change, 
certain internal structures of the interference pattern are constant because of the overall structure 
of the source. The correlation pattern observed with the Brown-Twiss stellar interferometer is 
a measure of this constant, internal structure. 
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Appendix A 


THE CONNECTION BETWEEN QUANTUM 
AND CLASSICAL DESCRIPTION 


Photons in the optical range are detected primarily by means of a detector that employs the 
photoelectric effect (e.g., a photoelectric tube). In fact, the detection of a photon in a light beam 
can be defined as the ejection of an electron from the photocathode of some suitable detector. In 
this context, the concepts of photon and classical electromagnetic field are related through a single 
consideration; that is, for a photocathode illuminated by a light beam, the instantaneous probability 
per unit time of emission of a photoelectron (detection of a photon) is proportional to the instan- 
taneous intensity of the electromagnetic wave falling on the photocathode. Thus, 

dP(l) = al(t)dt , (Al) 

where dP (1) is the probability for one, and only one, photon to be detected in a time dt; a is the 
proportionality coefficient. The probability for two photons to be detected is of second order in dt 
and is considered negligible. From probability theory, we know that for a finite interval of time 
between and t 2 , the probability for the detection of n photons is given by the Poisson distribution 


„ [ u (t i . t 2 )] n ex p [- u («= i . t 2 )] 

1 uO) " n ! 


(A2) 


where 


u 





al(t) dt 


Unfortunately, unless the intensity l(t) is a constant, the distribution P u (n) is not observable, 
for if one took statistics on successive intervals (t lt t 2 ), (t 2 , t 3 )---(t n ,t n+1 ), one would find 

that in general the value of u(t n , t n+1 ) would vary from interval to interval. Equation A2 would 
describe the distribution for only the collection of intervals that had a common value (e.g., 04 ) 
for the variable u (see Figure Al). 
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STATISTICS TAKEN ON ALL INTERVALS WILL NOT YIELD POISSON'S 
DISTRIBUTION DUE TO DIFFERENT VALUES OFu; a, , a 2 , a 3 , ETC. 


STATISTICS TAKEN ON THESE SELECTED INTERVALS WILL YIELD POISSON'S 
DISTRIBUTION SINCE u = a,ON ALL OF THESE INTERVALS. 


Figure A] — Sequence of time intervals having different values of u - 



l(t) dt. 


If the variable u varies from interval to interval such that the probability of any interval 
having a particular value of u is given by P(u) , the probability of detecting n photons in any 
interval selected at random will be given by the Poisson distribution averaged over all possible 
values of u; that is, 


P(n) 


P(u) 


u n exp (- u) du 
n ! 


(A3) 


It would be well at this point to ask what possible meaning can be given to the instantaneous 
intensity l(t) or its integral over an interval u(t n , t n + 1 ) since there seems to be no operational 
way to determine their values. The number of photons detected in any given interval has no unique 
relationship to the value of u for that interval. Therefore, I and u are calculations! devices used 
to determine the statistical properties of photoelectrons from a photocathode. If one calculates the 
probability per unit time for ejection of a photoelectron caused by a source of electromagnetic radia- 
tion some distance away and does so strictly from the theory of quantum electrodynamics, one 
finds that this probability is proportional to the absolute value squared of a vector quantity A. 
Furthermore, if the source is composed of a large number of elementary quantum systems and is 
some distance away, the quantity A may be calculated to a very good degree of approximation by 
using Maxwell’s equations and considering A to be the classical vector potential and hence I A 1 2 to 
be the classical intensity l(t) . Therefore, although the quantities i(u) and u(t n , t n + 1 ) are, 
strictly speaking, not observable (although their average values may be defined operationally), they 
are conceptually and calculationaily useful quantities. 
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Consider now a large number of nonoverlapping intervals of length T and the mean value and 
rms fluctuation of the number n of detected photons in each interval. Employing the definition of 
a mean value 


<f(n)> = £ f (")P(n) 

n=0 


(A4) 


and making use of the relationship 


co 

V -1 n! u n 

/ . (n - P )!irr ex p (~ u ) - uP 


(A5) 


we obtain 


00 /. CO 

(n) = J P(u) exp (-u)du 


j P(u)u 

J 0 


du - ^u) 


This may be taken to be the operational definition of (u), the average value of u 
The mean square fluctuation of n about its mean value will be given by 


(A6) 


(An 2 ) = <(n - ( n ) ) 2 ) 


= <" 2 > - <"> 2 

= (u 2 ) + <u> - (u ) 2 

= (n> + <(u - <u)) 2 ) 

= < n > + < Au2 > ■ 


(A7) 


The first term on the right-hand side of Equation A7 is the result that would be obtained if 
light were a stream of classical particles; the statistics of the photons then would be Poissonian. 
The second term (Au 2 ) , which is a measure of the fluctuations of the wave field caused by some 
sort of interference phenomena, shows that the deviation of photon statistics from Poissonian 
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because of their boson nature can be determined solely from an investigation of the associated 
classical electromagnetic field. In fact it can be shown (Reference 11) that the very existence 
of the classical electromagnetic field is intimately connected with the boson character of the 
photon. 

Consider now the situation where two different detectors A and B are situated at different 
points of space. For fixed values of u A and u B , the emission of electrons from the two cathodes 
is completely independent, giving 


P 

U A U B 


( n 


A’ 





(A8) 


and 


P ( n A- n B ) 





n B ) du A du B 


(A9) 


Equation A8 indicates the complete independence of the distributions of n A and n B for given 
values of u A and u B . However, if u A and u B are not distributed independently, i.e,, P(u A , u B j ? 

P(U A ) P(u B ), then P(n A , n B ) ? P(n A ) P(n B ) in general, and n A and n B are not distributed independently, 
In other words a correlation between u A and u B imposes a correlation between n A and n B . 

In the operation of the stellar interferometer of Brown and Twiss, the fluctuations in the 
"output" of two different detectors are multiplied together in a circuit called a correlator, and 
this product is time -averaged. The "output" of a photoelectric detector will be a current or 
voltage that is proportional to the number of photoelectrons ejected from its cathode during some 
interval of time T, which is essentially the resolution time of the detector. 

The averaged output of the correlator therefore is proportional to 


(( n A ' ( n A>) ( n B “ ( n B>)) " ( n A n B ) " ( n A ) < n B ) “ < U A U B>“ ( U a) ( U b) = ( A10 ) 

this output is directly expressible in terms of the fluctuations of the classical electromagnetic in- 
tensity. If the fluctuations of u A and u B are independent, then (u A u B }> = ^u A ^ , and there is 
no averaged output from the correlator. 
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Appendix B 

A DETAILED SOURCE MODEL 


Consider a model of a radiation source that is considerably more realistic than the one em- 
ployed in the section entitled "Refinements of the Model," Instead of constructing the source from 
a large number of oscillators with constant amplitudes and phases that change randomly over a 
period of time At , we shall construct the source from a collection of oscillators that turn on at 
random times with random starting phases and then after a time At turn off again. This model is 
suggested by the picture of a hot gas whose atoms are constantly being excited through collisions, 
resulting in the radiation of their extra energy during a short time At. The method to be used was 
employed by Rice (Reference 12) in his calculations of the "shot effect;" however, for a different 
treatment of the same problem the reader is referred to a paper by Janossy (Reference 13). 

Consider a collection of oscillators situated at a point s. If the i th oscillator at s turns on 
at t = o , it will produce a field at R given by 

A.(t) = e. exp(i0 i ) F(t - |SR|/c) , (Bl) 

where <£. is an overall phase angle, and e i is a unit polarization vector perpendicular to the line 
joining S and R. F(t) describes the characteristic output of an oscillator and has the properties 
F(t) = 0 for t < 0 and F(t) % o for T » At . Since the retardation factor |sr|/c will apply equally 
to all times at R, we shall drop it in subsequent calculations, remembering that any event at time 
t at s corresponds to an effect at R at time t + |sr|/c. 

The total field at R caused by many oscillators turning on at times t. with phase angle <f> i 
and polarization e. is given by 


A(t) = 


e i exp(i0.) F (t - t.) 


(B2) 


The field produced by N oscillators will be considered to be a random function of the 3n random 
variables <f> i9 , and t.. If the random variables have probability distribution P^), P(eJ , and 

P (t . ) , the average of a random function can be defined as the weighted integral of the function over 
all values of the random variables, 

<F(*i. «i. tl )> = IJ J P ^d^|p(*i)dS i |p(t l )dt l F(^,S i , t t ) . (B3) 
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The normalized autocorrelation function of the function F( t ) is defined as 


*(T) 


i; 


F* (t ) F(t + r)dt 


f* CO 

|F(t)l 2 

J-CD 


dt 


(B 4 ) 


where it is assumed that the integrals exist. It is easy to verify that $ (0) - l and $(-r) = <&*( T ) . 
Since F(t) is essentially zero everywhere except for 0 < t < At , $(r) is essentially zero every- 
where except for -At < r < At . 

Consider again the amplitude function, Equation B 2 . For a process that is stationary in time, 
the sum should be taken over an infinite number of oscillators over all times. Since it would be 
difficult to calculate averages with this quantity directly, we shall calculate with the subsidiary 
quantity A T (t ) , which is the amplitude produced by all of those oscillators with turn-on times t . 
such that -T/2 < t. < T/2. At the end of the calculation, we shall always take the limit T-<». 

It is assumed that the t . are independently and uniformly distributed with a probability per 
unit time 77, where 77T is the average number turning on in time T. The phase angles <£. are 
independently and uniformly distributed between 0 and 2 77 . The distribution of e. will be con- 
sidered later. 

The averages can be calculated in three steps; first, (f) tn is calculated for the case where 
exactly N oscillators turn on in time - T/2 < t < T /2 . We will then calculate 


<f) t = I> T ( N )<F) TN 

N=0 


where P X (N) is the Poisson distribution 


< :vDl [ 

N! 


exp ( 


V 1 ) 


Finally, 

( f} - xiS (f) t . 
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Now consider the intensity 


I(t) = A* (t) • A(t) 


(®i ‘ *i) ex P ~ (B5) 


(*( t ))nt 




N 

X (ej • 6j) exp(i^>j - i^>i) F*(t - tj) F(t - tj) 

_ iTFi 


(B6) 


Only those terms in the sum where i = j give a non- zero contribution after the phase angle 
averaging; therefore (recalling e. • e ; = l), 


<Kt)> NT = 



(B7) 


Since the time interval T can be as large as desired, as long as t is not too close to 
-T/2 or T/2 , we may replace the integral over | F | 2 by the infinite integral to obtain 


<KO> nt 



I F( t ) | 2 dt 


(B8) 


Also, 


<I(t)> T = 



(B9) 


The passage to the limit T - “ is now trivial, and thus 


<KO> 



I F( t ) | 2 dt 


(BIO) 
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The average intensity is just the intensity produced by a single oscillator multiplied by the 
average number of oscillators turned on per unit time-a rather intuitive result. The result is 
independent of the state of polarization, i.e., independent of I s (e i ) . 


Consider now the autocorrelation function of the intensity (i//(r)) defined as \p(r ) = 

<l(t) I (t + r)>: 


Wr)\ 


n r 2n d <j>. r f T/2 dt . 

n 1 L ^ 


N 

Y, (e; • e.) (e k • exp i(^ j - + 4>l ~ <£ k ) 

_j . j >k,£=l 


! (t _t i) F (t + T " h) F (t + T_ h) 


Again phase angle averaging gives a zero result unless i = j and k = l , or i = l and j = k, 
or i = j = k = l . This gives 


OK -OX 


N r f T/ 2 rT 

Y T lF(t - ti ) | 2 dt 4 Y J 

ife=i L J ' T/2 JL 1 


|F(t + r-t k )| 2 dt k 


Y P (®k)(®i • £ k ) 2 d ®i d «k 


rr/2 rT/2 

T F*(t- ti )F(t + r- ti ) dtj T F *( 

» — T / 9 J “T/2 


Fit ~ tk) dt i 


N rT/2 

r I iF(t - 1.)\ 2 iF( t + r - t i )i a dt i 

L J J - t /2 


ii ii ■■ll■llll■llllln■lll i n i « n m in ii i hi mu \m hi mi ll■l■ll«ll[ I II 1 1 II in 


Vl_ 


Converting the integrals to infinite integrals under the assumption that both t and t + r are 
far removed from -T/2 and T/2, gives 


[«r)] NT 


N(N - 1) 
T 2 



I F(t ) | 2 dt 


N(N - 1) 
T 2 


<(•. 



F*(t)F(t + r) dt 


2 


N 

+ T 



I F( t ) | 2 [ F( t + r) | 2 dt 


(B13) 


where 


<(e, ' e k ) 2 > = ffp(e ; ) P(« 


<)(« 


i , ) 2 de . de 

k J i r 


Averaging over N and taking the limit T-°o gives 


r) = V 2 


f* 00 

|F(t)| : 

J — 00 


dt 


V — o 


F*( t ) F(t + r) dt 


<(*, ' *•.)’> 


/* CD 

|F( t ) i 2 I F( t + r) | 2 dt 

» -00 


(B14) 


Turning now to the quantity ([e i ■ e k ) 2 ^> we define two states of polarization, completely 
polarized by P(e.) = s(ej - e') and completely unpolarized by P (e .) = (27t) _1 . The quantity 
<(ej • e k ) 2 ) = i for completely polarized and ((e. ■ e k ) 2 ^> = 1/2 for completely unpolarized. 
Equation B14 can be simplified to be 


Wt) 




|F( t ) I 2 | F( t + r) | 2 dt 


(B15) 
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r 


where y - 1 for unpolarized and = 2 for polarized. 

In the limit of a very large tj (equivalent to a large number of oscillators being on at any given 
time), we may neglect the term linear in -n compared to (i) 2 , which is quadratic in r), giving 

iKt) = <l ) 2 [l l$(r)| 2 ] , ( B16 ) 

a familiar result (Reference 14). For r = 0, <D(t) = 1; this gives 

\/-(0) = <I 2 (t)> = (l+^)<I(t)> 2 . (B17) 

our desired result. 
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